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We study systems analogous to binary black holes with spin in order to gain some insight into the 
origin and nature of "bobbing" motion and "kicks" that occur in this system. Our basic tool is a 
general formalism for describing the motion of extended test bodies in an external electromagnetic 
field in curved spacetime and possibly subject to other forces. We first show that bobbing of exactly 
the type as observed in numerical simulations of the binary black hole system occurs in a simple 
system consisting of two spinning balls connected by an elastic band in flat spacetime. This bobbing 
may be understood as arising from the difference between a spinning body's "lab frame centroid" 
and its true center of mass, and is purely "kinematical" in the sense that it will appear regardless of 
the forces holding two spinning bodies in orbit. Next, we develop precise rules for relating the motion 
of charged bodies in a stationary external electromagnetic field in flat spacetime with the motion 
of bodies in a weakly curved stationary spacetime. We then consider the system consisting of two 
orbiting charges with magnetic dipole moment and spin at a level of approximation corresponding 
to 1.5 post-Newtonian order. Here we find that considerable amounts of momentum are exchanged 
between the bodies and the electromagnetic field; however, the bodies store this momentum entirely 
as "hidden" mechanical momentum, so that the interchange does not give rise to any net bobbing. 
The net bobbing that does occur is due solely to the kinematical spin effect, and we therefore argue 
that the net bobbing of the electromagnetic binary is not associated with possible kicks. We believe 
that this conclusion holds in the gravitational case as well. 



I. INTRODUCTION 

The successful simulation of binary black holes with 
spin [1-4] has presented two intriguing surprises: If the 
spin lies in the orbital plane, an up-and-down "bobbing" 
of the entire orbital plane is seen to occur in phase with 
the orbit, and when the black holes merge, the final black 
hole can acquire a "kick" velocity exceeding those occur- 
ring in the non-spinning case by an order of magnitude. 
In the kick-maximizing case of circular orbit with an- 
tiparallel spins lying in the orbital plane, the direction 
(up or down) and speed of the kick depend on the or- 
bital phase at merger. Two natural questions are: What 
causes the bobbing motion? Can the kick be viewed as 
a post-merger continuation of the bobbing? 

In this paper we shall address these questions by study- 
ing simple analogous systems: analytical toy models of 
the spinning black hole binary system that display some 
of its salient features. We first consider two spinning bod- 
ies connected by an elastic band in flat spacetime, a sys- 
tem we refer to as "tetherballs." We find that this system 
displays bobbing motion very similar to that observed in 
numerical simulations of binary black holes. The bob- 
bing is identified as arising from a purely "kinematical" 
spin effect that gives rise to bobbing independently of 
the details of the particular type of force holding the two 
bodies in orbit. This effect is analogous to the Thomas 
precession, which gives rise to a precession independently 
of the details of the particular type of force accelerating 
a spinning body. 

The analog system of primary interest to us is an elec- 
tromagnetic binary, consisting of two charged magnetic 
dipoles with spin in fiat spacetime. Here we find that, 
in addition to the "kinematical bobbing" associated with 



their spin, the individual bodies also undergo "dynamical 
bobbing" associated with their magnetic dipole moment. 
However, at the level of approximation corresponding 
to 1.5 post-Newtonian — which we refer to as "1.5 post 
Coulombic" — we shall show that the mass- weighted dy- 
namical bobbing of the two bodies is equal and opposite, 
so that the net bobbing (i.e., the sum of the mass times 
the bobbing velocity for the two bodies) is purely kine- 
matical, just as in the case of tetherballs. Nevertheless, 
we will also see that a considerable amount of momentum 
is exchanged between the bodies and the electromagnetic 
field during each orbit. This momentum exchange does 
not contribute to the net bobbing because the net mo- 
mentum acquired by the bodies from the electromagnetic 
field goes into the "hidden mechanical momentum" [5, 6] 
of the bodies. 

If an electromagnetic analog of a merger event were 
to occur, a net kick could be experienced by the system 
only if momentum is carried away to infinity by the elec- 
tromagnetic field. We cannot analyze the momentum 
radiated to infinity — momentum cannot be radiated to 
infinity at the 1.5 post-Coulombic approximation — but, 
in a more exact treatment, we would expect the radiated 
momentum to be related to the total momentum stored in 
the electromagnetic field at or near merger. However, the 
momentum stored in the electromagnetic field is entirely 
unrelated to the net bobbing of the bodies. In partic- 
ular, the momentum stored in the electromagnetic field 
depends upon the magnetic dipole moment of the bodies 
but not their spin, whereas the net bobbing depends upon 
their spin but not their magnetic dipole moment. Thus, 
in the electromagnetic case, we believe that the bobbing 
phenomenon is entirely unrelated to possible kicks. 

In this paper, we will derive precise rules for relating 
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the equations of motion of charged bodies in a stationary 
external electromagnetic field in flat spacetime with the 
motion of bodies in a weakly curved stationary space- 
time. These rules will enable us to directly translate our 
electromagnetic results on bobbing motion into results 
for the corresponding gravitational case — i.e., a binary 
system of spinning masses — at the 1.5 post-Newtonian 
approximation. This gravitational binary system has 
been directly analyzed in the 1.5 post- Newtonian ap- 
proximation in [8], and our results on bobbing in the 
gravitational case agree with theirs. However, it is much 
more difficult to unambiguously interpret the results in 
the gravitational case as compared with the electromag- 
netic case for the following three reasons: (i) In the grav- 
itational case, the kinematical bobbing terms take the 
exactly same mathematical form as the dynamical bob- 
bing terms, so it is much less clear how to make a clean 
distinction between effects of a kinematical origin and ef- 
fects of a dynamical origin, (ii) The description of the 
motion of bodies in the gravitational case is gauge depen- 
dent, as are the "gravitational field strengths." (iii) In 
general relativity, only the total momentum of the entire 
system is well defined. Any allocation of a portion of the 
momentum of the system to the gravitational field and 
a portion to the bodies requires both a choice of pseu- 
dotensor (such as the Landau-Lifshitz pscudotensor used 
in [8] ) and a choice of gauge, and it is unclear the extent 
to which this allocation may depend on these choices. 

Nevertheless, the very close analogy established in 
this paper between the electromagnetic and gravitational 
cases leads us to believe that our interpretation of the 
nature of bobbing and kicks in the electromagnetic case 
should extend to the gravitational case. We therefore 
propose the following answers to the questions posed in 
the first paragraph above: The bobbing seen in simula- 
tions of black hole mergers is a purely kinematical spin 
effect that is in no way special to the gravitational in- 
teraction. The kick has an entirely different cause, and 
should not be viewed as a post-merger continuation of 
the bobbing phenomena 1 . 

The nature of the effects we wish to analyze in this 
paper require us to be very careful and precise in our 
definitions of quantities such as the linear momentum, 
spin, and center of mass position of a body, which will 
require us to develop a substantial amount of machin- 
ery in the next section. It may therefore be worthwhile 



As previously mentioned, in the electromagnetic case the net 
bobbing depends upon spin of the bodies, whereas the momen- 
tum stored in the electromagnetic field (which we have associated 
with kicks) depends upon magnetic dipole moment of the bod- 
ies. Under the translation rules between electromagnetism and 
gravity established in section IV, these distinct terms in electro- 
magnetism translate into terms of exactly the same functional 
form in gravity. Therefore, it should not be surprising if in the 
gravitational case, the magnitude of the kick velocity and its de- 
pendence on the phase of the orbit are very similar to those of 
the net bobbing velocity. 



to give here a simple explanation of the nature and ori- 
gin of the kinematical bobbing effect in flat spacetime. 
Consider a body described by stress-energy tensor T M „ in 
global inertial coordinates (t,x l ) that is acted upon by 
an external 4- force density / M , so that T M „ satisfies 
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Integrating this equation, we obtain 



~ / T 0u d 3 x = / f u d 3 x. 



(1) 



(2) 



Multiplying the v = component of Eq. (1) by x" 1 , in- 
tegrating over space, integrating by parts, and using (2), 
we obtain 



dt 



= / T 



+ f(x*-2i)d*X. 



(3) 



where £ = J T 00 d 3 a; is the total energy of the body as 
measured in the lab frame, and we have defined the "lab 
frame centroid" 



~ / T"".r'd '.r. 



(4) 



Differentiating (3) with respect to time and using Eq. (2) 
again, we obtain 



£ 



iM 
dt 2 



fd 3 x 



d 



/V 
dt 



zl)d 3 x 
fd 3 x. 



(5) 



We now restrict consideration to non-relativistic mo- 
tion and drop terms of order |d,?L/di| 2 . We may then 
replace £ by the rest mass 2 , to, of the body. If the in- 
tegrated 4-force density is nearly orthogonal to the 4- 
vclocity of the lab frame centroid (as must occur if the 
rest mass of the body is approximately constant), then 
the last term on the right side of (5) will be of order 
|dzL/dt| 2 and may be dropped. Furthermore, if the ex- 
ternal force is nearly centered about the lab frame cen- 
troid, the second term on the right side of (5) also will 
be negligible 3 . Dropping this term as well, we obtain 



dt 2 



fd 3 x, 



(6) 



2 The rest mass of a body is given by m = (— p A 'p M ) 1/ ' 2 , where 
p^ is the 4-momentum of the body, defined (in the absence of an 
electromagnetic field) by an integral of T^ v over a spacelike slice. 
However, even in Minkowski spacetime, if an external force acts 
on a body, then p^ will be slice dependent. In particular, p^ and 
m depend on the choice of frame and we must distinguish be- 
tween, e.g., the rest mass as determined in the lab frame and 
the rest mass mcM as determined in the center of mass frame. 
These differences will be negligible for the present calculation, 
and we will ignore them. Similar remarks apply to spin. We will 
treat these issues with great care in section II. 

s If this term is not negligible, it may contribute to additional 
"dynamical bobbing" of the body. This occurs in the electro- 
magnetic case. 
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FIG. 1. An Illustration of kinematical bobbing. The lab 
frame centroid of a body (solid line) is confined to the x x -x 2 
plane. However, since the body has a (constant) spin vector 
lying in the plane, its center of mass (solid dot) is displaced 
in the x 3 direction by a velocity dependent factor (equation 
(7)), giving rise to bobbing. The body is shown at two phases 
in its orbit. In the foreground, the bottom portion of the 
body contains extra lab frame kinetic energy (illustrated with 
shading), skewing the lab frame centroid downwards so that 
the center of mass must have bobbed up. One half of an orbit 
later (shown in the background), the situation is reversed and 
the center of mass has bobbed down. If two bodies with 
anti-aligned spins were considered, a net bobbing would be 
observed. 

which can be recognized as a version of F — ma. How- 
ever, the "a" that appears in this equation is the accelera- 
tion of the lab frame centroid. The lab frame centroid z^ 
measures the "center of energy" of a body as determined 
by a lab frame observer. However, for a spinning, mov- 
ing body, this center of energy will be skewed towards the 
side of the (rotating) body that moves with the body's 
bulk motion, on account of the extra lab frame kinetic 
energy present there. As obtained in equation (60) be- 
low, the true center of mass z l CM is related to the lab 
frame centroid approximately by 

mz C M = mzl + Sx ^+0(|^| 2 ). (7) 

If J f l d 3 x lies in a plane at all times, and if dz^/dt ini- 
tially lies in this plane, then by (6) the lab frame centroid 
is confined to the plane. However, if the spin of the body 
has a component in the plane, then by (7) the center of 
mass will be displaced out of the plane. This accounts 
for the "kinematical bobbing" (see figure 1). 

The remainder of this paper is organized as follows. In 
section II we present a general formalism for describing 
the motion of extended test bodies in curved spacetime. 
We explicitly include the possibility of an external elec- 
tromagnetic field, and allow for the possibility of addi- 
tional external forces. The first subsection defines linear 



4 A derivation of equation (7) is also given in [8] (see also [25, 26]). 
However, the notion of center of mass used in [8] corresponds to 
the "hybrid" definition of center of mass introduced in subsection 
II B and does not agree with the notion of center of mass used 
in the present paper. 



and angular momenta of a body. The second subsection 
deals with the association of a center of mass worldline to 
a body. The third subsection deals with the evaluation 
of electromagnetic and gravitational forces and torques 
up to quadrupole order. Section II is almost entirely re- 
view, although some aspects of the presentation are new 
and there is some new material on, e.g., the relationship 
between the true center of mass and the "lab frame" and 
"hybrid" centroids. 

In section III we apply this formalism to derive bobbing 
effects in simple example systems. First we consider the 
case of tetherballs and derive the kinematical bobbing 
effect discussed above. We then consider the example of 
a magnetic dipole orbiting a line electric charge in flat 
spacetime, and a spinning body orbiting a line mass in 
general relativity. In these cases the kinematical bobbing 
is supplemented by dynamical bobbing. 

In section IV we give a careful treatment of the anal- 
ogy between the equations of motion of a body mov- 
ing in the electromagnetic field of a stationary source 
in flat spacetime and the motion of a body in a weakly 
curved, stationary spacetime in general relativity. We 
extend the results of [7] to include quadrupole terms and 
terms quadratic in the velocity of the body. We provide 
precise rules for translating the electromagnetic evolu- 
tion equations for center of mass motion and spin into 
the corresponding gravitational equations at this order 
of approximation. 

Finally, in section V we consider the electromagnetic 
analog of a black hole binary, consisting of two orbit- 
ing charged magnetic dipoles with spin. We assume that 
each body moves in the electromagnetic field of the other 
according to the laws of test body motion derived in sec- 
tion II. (These laws were also obtained in [9, 10], where 
self- force effects were considered in addition.) We ob- 
tain the spin and magnetic dipole terms appearing in 
the equations of motion at "1.5 post-Coulombic" order 
(corresponding to 1.5 post-Newtonian order in the grav- 
itational case). Using the translation rules developed in 
section IV we also obtain the corresponding results for 
gravity, which agree with [8]. We analyze the flow of 
momentum between the bodies and the electromagnetic 
field, and discuss what can be learned from this analysis 
about kicks. 

Our notation and conventions on indices are as follows. 
Latin indices from early in the alphabet (a, 6, c, . . . ) are 
used to denote abstract tensor indices. Greek indices are 
used to denote components of these tensors (most com- 
monly in global inertial coordinates in Minkowski space- 
time). Latin indices from the middle of the alphabet 
(i, j, k, . . .) denote spatial components, and a "0" will be 
used for time components. Sign conventions are those of 
[11]. 
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II. MECHANICS OF EXTENDED TEST BODIES 

Consider a body with stress-energy tensor T ab — T^ ab ^ 
and electromagnetic charge-current density J a moving 
in a curved spacetime (M,g ab ). Both T ab and J a are as- 
sumed to have support within a spatially-compact world- 
tube W. We also allow for an electromagnetic field F ab , 
which is assumed to satisfy Maxwell's equations 

V [a F bc] = (8a) 
V b F ab = 4mJ a (8b) 

in a neighborhood of W . It follows from this that the 
current density is conserved: V a J a = 0. 

The total stress-energy tensor T^ b t of the system is also 
assumed to be conserved. (This would follow from Ein- 
stein's equation, but we do not impose Einstein's equa- 
tion in the present section.) T^ b t is the sum of the body 
stress-energy T ab , the stress-energy of the electromag- 
netic field 

^em = ^(F ac F b c - \9 ah F cd F^), (9) 

and the stress-energy T£ b t associated with any other fields 
or bodies that may be present. Conservation of this total 
yields 

V fc T ab = F ab J b + f a , (10) 

where F ab J b = — V^T^m is the ordinary Lorentz force 
density and f a = —V b T£ b t is the "external force density." 
Except in Sect. Ill A below, this paper will mainly be 
concerned with the case f a = 0. 

We now proceed to review some basic results on the 
mechanics of compact bodies satisfying the equations just 
described. Linear and angular momenta are defined in 
Sect. II A, and various aspects of their evolution are dis- 
cussed there. The notion of a center of mass worldline is 
then introduced in Sect. II B. For stationary spacetimcs, 
we also define the notion of a lab frame centroid and 
compare it with the center of mass. Multipole approxi- 
mations for the force and torque are presented in Sect. 
II C up to quadrupole order. The formalism described 
here is mainly due to Dixon [12-14] (see also [15-19]), 
although many elements of our presentation arose more 
recently in [20]. 

In order to be mathematically well-defined, the defini- 
tion of momentum given in Sect. II A requires only that 
W lie in a suitable normal neighborhood of the timelike 
worldline used in the construction of certain "general- 
ized Killing fields." This imposes relatively weak restric- 
tions on the size of the object and the curvature within 
it. It does not, however, restrict one to the test body 
regime. Electromagnetic self-fields may be large, and, if 
Einstein's equation is imposed, the body can contribute 
significantly to the spacetime curvature. Some relatively 
minor additional restrictions are required in Sect. II B 
in order to ensure the existence of a unique center of 



mass worldline [21], but the constructions and results 
presented there also do not require that any self-fields be 
negligible. 

This contrasts with the multipole approximations 
given in Sect. II C, which are valid in the form presented 
only for sufficiently small test bodies. Despite this, it can 
be shown that the definition of momentum given in (24) 
below can be appropriately modified so that essentially 
identical approximations remain correct even in the pres- 
ence of large self-fields [9, 22, 23]. Indeed, using a modi- 
fied definition of the quadrupole and higher moments of 
T ab (but not the moments of J a ), the laws of motion are 
changed only in that each instance of the external field is 
replaced by that external field plus the so-called "regu- 
lar" Detweiler- Whiting component of the body's self-field 
[24]. A short review of the full theory may be found in 
[23] for the special case of self-interaction with a scalar 
field in a fixed curved spacetime. This generality will not 
be needed here, so we shall adopt definitions equivalent 
to Dixon's for simplicity. Consequently, the applicabil- 
ity of the results of this section should be considered as 
restricted to the test body regime. 

A. Momentum 

Standard definitions for a body's linear and angular 
momenta in flat spacetime may be found in textbooks 
(see, e.g., [25, 26]). In global inertial coordinates a;^, 
they are given by 

pt* = [ T^dS u (11a) 

S^(z) = S [ ^\z) = 2 [ (x- z)^T^ x dS x , (lib) 

where is an arbitrarily chosen origin in spacetime and 
£ any Cauchy hypersurface. It is common to also intro- 
duce a spin angular momentum vector 

S, = -\e^ p [ 7 £=) S^z). (12) 

Assuming that F ab J b + f a = 0, it immediately follows 
from (10) that p^, S^ v , and do not depend on £ for 
systems with finite radius. p M and are independent of 
z M as well. 

It is useful to rewrite (11) without any reference to 
particular coordinate systems. The resulting procedure 
combines both the linear and angular momenta into a 
single expression, and makes explicit that they are both 
associated with Poincare symmetries. Given any Killing 
field f, let 

V e = f T a b £ b dS a . (13) 

■pj may be viewed as a linear map from the space of 
Killing fields into R. Since any Killing field is uniquely 
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characterized by its value and that of its first derivative 
at an arbitrarily chosen point z, we may express as 
a linear map on the space of "Killing data" £ a (z) and 
V a £ fc (z). It therefore has the form 

^=/(z)^(z) + ^ ttfc (z)V fl 6W (14) 

for some tensors p a (z) and S ab = S^ ab \z) defined at z. 
These are exactly the linear and angular momenta given 
by the more familiar definition (11). Note that for fixed 
£ a , is independent of both z and £ when F a bJ b + f a = 
0. It therefore generates a class of conserved scalars in 
this case. 

Eqs. (13) and (14) are trivially extended to apply in 
any maximally symmetric spacetime. As we shall now 
show, a suitable generalization of the notion of a Killing 
field allows them to be applied in generic spacetimes [20] . 



1. Generalized Killing fields 

The appropriate notion of a generalized Killing field 
will depend upon the choice of a smooth timelike world- 
line Z = {z(s)|Vs € M.} and a unit future-directed time- 
like vector field n a (s) G T z r a \M defined along Z. Note 
that s is not required to be proper time and n a need not 
lie tangent to Z. A "center of mass frame" will be de- 
fined in Sect. II B below that uniquely fixes Z and n a (s). 
For the constructions and definitions of this subsection, 
however, they may be chosen arbitrarily, subject only to 
the body being contained in the normal neighborhood of 
Z defined below. 

A specific generalized Killing field may be fixed by 
choosing a time s together with tensors A a (s) and B a b — 
B[ab](s) at z(s). We use the Killing transport equations 



D 

ds 
D 

d^ 



A a (s) - z b { s )B ba {s) = (15) 
B ab {s) + R abc d {z{s))z c {s)A d {s) = (16) 



to uniquely extend these tensors to all of Z, where 
D 



d.s 



i a V a . 



(17) 



Note that the skew symmetry of B a b is preserved by this 
prescription . 

Now consider all pairs {z, v a ), where z € Z and v a € 
T Z M is orthogonal to n a . This forms a subset T±Z of the 
tangent bundle TM. For any element of T±Z, one may 
associate an affinely-parameterized geodesic y(u) whose 
initial point is y(0) = z, and whose initial tangent is v a . 



5 It is possible to use initial data for which B a t is not symmetric. 
The vector field that eventually results generalizes a homothety 
or other non-Killing affine collineation [20]. 



As long as these geodesies can be extended sufficiently 
far, the map (z,v a ) — ¥ y(l) will be a smooth function 
from T±Z to M. Its Jacobian is clearly invertible at 
every point (z,0), so it follows from the inverse function 
theorem that the given map defines a diffeomorphism on 
some neighborhood U C M oi Z. This will be the region 
in which the generalized Killing fields (GKFs) are to be 
defined. 

We now define the GKF £, a (x) in U associated with 
our choice of initial data A a {s) and B a b(s). Let x e U. 
The diffeomorphism described in the previous paragraph 
may be used to uniquely associate x with an element 
(z(t x ), v a {r x )) of T±Z. Use this pair to construct a 
geodesic y(u) as before. The GKF may then be computed 
along y(u) by solving the Jacobi (or geodesic deviation) 
equation 



D 2 e a 

du 2 



Rabc d y b y c U = 



with initial data 



£,a{z{r x )) = A a (r x ) 



du 



v b B ba {r x ) 



(18) 

(19) 
(20) 



Here, y a denotes the tangent to y(u). The given equa- 
tions uniquely define £ a (a0 throughout U once A a and 
B ab are given at any one point on Z. 

It will often be useful to introduce the spacelike hy- 
persurface S(s) C U as the image of (z(s) 7 v a ) as v a is 
varied while s is held fixed (i.e. the union of all geodesies 
mU passing through z(s) and orthogonal to n a (s) at that 
point). Restricted to any one of these hypersurfaces, £ a 
is an ordinary Jacobi field. All solutions to the Jacobi 
equation can be written down explicitly in terms of the 
first and second derivatives of Synge's world function as 
well as the initial data A a and B a b [12]. The first of these 
two tensors is easily interpreted using (19): It specifies 
£ a on Z. Combining (15) and (20) shows that B a b serves 
to fix 



V a £ & (z(s)) =B ab {s) 



(21) 



for all s. 

One way to motivate the given definition of the GKFs 
is to rewrite (18) as a pair of first order differential equa- 
tions. The result is very similar to the Killing transport 
equations: An analog of (15) carries over exactly, while 
(16) is modified only by an overall contraction with y b . 
It is necessary to modify the Killing transport equations 
in this way because inconsistencies might arise if one at- 
tempts to propagate derivatives transverse to the direc- 
tion of transport along all spatial geodesies emanating 
from a single point. Further details and alternative mo- 
tivations for our definition of generalized Killing fields 
may be found in [20]. 

The space V(Z, n) of all possible GKFs associated with 
a particular choice of Z and n a (s) is obtained by vary- 
ing the initial data A a {s) and B a b(s). A general GKF is 
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clearly linear in this data, so V(Z, n) is a vector space. In 
four spacetime dimensions, dim V(Z, n) — 10. Also note 
that specifying a GKF and its first derivative at any one 
point on Z fixes it throughout U. Since any true Killing 
field satisfies both the Killing transport equations along 
any curve and the geodesic deviation equation along any 
geodesic, it is clear that true Killing fields are also gen- 
eralized Killing fields. In a maximally symmetric space- 
time, V(Z, n) is exactly the space of all Killing fields 
(and therefore independent of Z or n a (s) in such a case). 
Finally, we note that for any £ a £ V(Z,n) and z € Z, 
one may show that 



££9ab{z) = V a C^g bc (z) = 0. 



(22) 



From this, it follows by the same derivation as for true 
Killing fields (see, e.g., [11]) that 



V„V 6 £ c (aO = £, d {z)R dab c (z). 



(23) 



for any z <E Z. 



2. Definition of momentum 

We may now define the momentum of an essentially 
arbitrary compact object moving in a general spacetime. 
Choose a timelike worldline Z and a vector field n a on Z 
as above. We require only that the body's worldtube W 
be contained within the neighborhood hi in which gen- 
eralized Killing fields are defined (at least for the time 
interval under consideration). The definition (13) for the 
momentum could now be carried over by simply replac- 
ing the arbitrary Killing field appearing there with an 
arbitrary GKF. It is, however, useful to add an addi- 
tional term in the presence of an electromagnetic field 
[12, 14, 27]. For the remainder of this paper, we define 
■Pj at time s by 



V£s)= / dS a T\(x)e(x) + J a (x) 
Ms) L 

x / duti-VV.^G/O^Mi/') 

Jo 



(24) 



The £ a appearing here is an arbitrary generalized Killing 
field. S(s) is the aforementioned spacelike hypersur- 
face generated by geodesies orthogonal to n a (s) at z(s). 
The u-integral is taken along the affinely-parameterized 
geodesic y'(u) satisfying y'(0) — z(s) and y'(l) — x. Also 
note that the biscalar a is Synge's world function. Its 
derivative a a >{y',z) = V a 'a(y' , z) is well-known to lie 
tangent to y'(u). 

The additional electromagnetic term appearing in (24) 
as compared with (13) is useful for constructing con- 
servation laws in symmetric systems: See [12] or Eq. 
(29) below. It also arises from the general theory of 
multipole moments describing J a and T ab . At least if 
f a = 0, there is a sense in which evolution equations for 



V% as defined above are the only ones implied by stress- 
energy conservation [14]. The associated linear and an- 
gular momenta — see (25) below — are also the canonical 
choices in a Lagrangian formalism describing the motion 
of charged particles with respect to some reference world- 
line (such as Z) [27]. However, in a multipole expansion, 
the dominant contribution of this additional electromag- 
netic term will enter via the electric dipole moment of 
the body. Since we will be interested below only in the 
case of bodies with vanishing electric dipole moment, this 
term will play essentially no role in the remainder of this 
paper. Nevertheless, we include it here for completeness. 

Since 'Pj(s) is linear in £ a , and each GKF is uniquely 
determined by its value and that of its derivative at z(s), 
we may write 



(25) 



for some p a and S ab — S>- ab '. It is not difficult to see 
that these tensors depend only on the single point z(s) 
and vector n°(s); the extension of Z and n a away from 
z(s) does not affect any GKF restricted to S(s), and is 
therefore irrelevant for the definitions of p a and S ab at s. 
Although presented differently, the definitions (24) and 
(25) for p a (z,n) and S ab (z,n) are equivalent to the ones 
originally given by Dixon [12, 14]. 

We now derive equations for the time evolution of the 
momenta. For this purpose, it is convenient to rewrite 
(24) in terms of a vector potential A a satisfying F ab = 
2V[q^4m. Integrating by parts, we obtain 



£(s) 



dS a \T\(x)e(x) + J a (x) 



(M b )\l 



duu 1 a b (y,z)C^A b (y) 



(26) 



Differentiating this using (10) yields 

'1 

di> c w" 



_d_ 

d^ 



(V^+qAaC) = f dS c 



T ab C i9ah + U a 



J a U s A a - V Q J du u- x a b C^A b 



(27) 



where w a is the time evolution vector field for the hyper- 
surfaces E(s). The combination A a £ a on the left-hand 
side is understood to be evaluated at z(s), and q is the 
total electric charge 



1 



J a dS a . 



(28) 



One important special case of (27) arises in the pres- 
ence of an exact Killing vector ( a that also preserves the 
electromagnetic field, i.e., C^g a b = £-QF ab = 0. It is then 
possible to choose an A a such that £(A a = 0. In such a 
gauge, it is clear that the right-hand side of (27) vanishes 
when f a — 0. This means that 



V c ( S ) + q(CA a )\ z(s) 



(29) 
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is a constant of the motion. In the absence of an electro- 
magnetic field (or if q = 0), Vc, is conserved by itself. 

Another interesting application of the constancy of 
(29) is to the case of a charged particle in flat space- 
time with an electromagnetic field satisfying V a F bc = 0. 
If £ a is a purely translational Killing field (so V a C fc = 0), 
one finds that 

^ = e^ = q eF ab z b . oo) 

ds ds 

This has the form of the Lorentz force law, and is valid 
without any restriction on the size of the body. It should, 
however, be noted that at this stage, Z is an arbitrary 
worldline. Although the momentum Vq evolves along Z 
by the same formula as the momentum of a point charge 
obeying the Lorentz force law, the center of mass world- 
line of the body (defined in the next subsection) does not 
necessarily coincide with a Lorentz force trajectory. 

Returning to the general case without symmetries, we 
may obtain general evolution equations for p a and S ab by 
substituting (25) into (27) while using (23). This gives 

-^=--R a bcd z b S cd + F a (31a) 
ds 2 

-4— - 22>i b J + N ab , (31b) 
ds 

where F a and N ab — are tensors at z(s) obtained 

by using (27) to express dV^ /ds in terms of the data for 
the generalized Killing field £ a at z(s) via 

^ = F a £ a + ^N ab V a ^ b . (32) 

We refer to F a and N ab , respectively, as the net force 
and torque acting on the body. Knowledge of dV$ / ds for 
all generalized Killing fields £ a is completely equivalent 
to knowledge of F a and N ab . 

When the force and torque vanish, (31) reduces pre- 
cisely to the equations derived by Papapetrou to describe 
the motion of a spinning (but otherwise "structureless") 
test particle. The Papapetrou terms in our equations 
arise from how the data for a generalized Killing field 
vary with z(s), and are thus "universal" in the sense that 
they do not depend on the details of the definition of 
V$ . From (27), it is clear that both the force and torque 
vanish entirely in maximally symmetric spacetimes with 
F ab = fa = 0. In this case, the Papapetrou equations are 
exact for arbitrary extended test bodies. More generally, 
approximations for the force and torque will be given in 
Sect. IIC below. In a multipole series, the lowest order 
expressions for the gravitational force or torque both in- 
volve the quadrupole moment of the body's stress-energy 
tensor. By contrast, the electromagnetic torque includes 
dipole and higher moments of J a . The electromagnetic 
force, of course, starts at monopole (Lorentz) order. 



B. Center of mass and motion 

As it stands, the formalism just described leaves Z and 
n a (s) arbitrary. We now fix these objects in such a way 
that Z can be interpreted as the body's center of mass 
worldline. 

Under mild assumptions, it may be shown that there 
exists a unique future-directed unit timclike vector 
riQ M (z) associated with every point z such that 

p a {z,n C M(z)) = m C M(z)nc M (z) (33) 

for some "mass" mcM( z ) > [21]. Furthermore, there 
exists a unique timelike worldline formed out of all points 
zcm satisfying 

P a {zcM,n C M(zcM))Sab(zcM,n C M(zcM)) = 0. (34) 

The resulting worldline and vector field define what we 
shall call the center of mass frame for the given body 5 . 
Roughly speaking, (33) defines the preferred hypersur- 
faces of "constant time" to be those attached to an in- 
stantaneously zero-momentum "rest frame." Similarly, 
(34) essentially implies that the body's mass dipole mo- 
ment about zcm vanishes in this same frame. Unless 
otherwise stated, the given definitions for the reference 
worldline and foliation will be assumed throughout the 
remainder of the paper. "CM" subscripts will be dropped 
for brevity when no confusion can arise. 

An equation of motion for the center of mass worldline 
z(s) may be obtained by differentiating (34) with respect 
to s and using methods introduced in [18]. In doing so, 
it is convenient to parameterize Z such that n a z a = — 1. 
Note that this implies that s is not a proper time along 
Z unless n a is tangent to Z. The result is that 

mi a =p a - N a b n b 

S ab [mF b + (p c - N c d n d )(qF bc - \S d f R bcdf )] 

m*-\S b c{qF bc -\S d lR bcdf ) ' 1 ' 

The m appearing here was defined by (33), so m = 
\J—VaV a - We have also introduced a "reduced force" 
F a = F a - qF a b z b . 

The fact that the momentum p a fails, in general, to 
equal mz a is fundamental to the results of this paper. 
When p a mz a , the body may be said to possess "hid- 
den momentum." This concept has arisen in various con- 
texts before (see, e.g., [5, 6, 28-30]), although it does not 
appear to have previously been discussed with this gen- 
erality. Note that the existence of hidden momentum can 



6 An alternative possibility that may be found in the literature 
sets n a <x z a and z a S a t, = 0. This is not, however, an acceptable 
choice. In general, the resulting "center" of even a free mass in 
flat spacetime would be found to satisfy an equation of motion 
that is third order in time. Its solutions include unphysical helical 
trajectories as well as the expected straight-line motions [31] 



introduce some confusion regarding what is meant by the 
"force" on a body. As is clear from (31a), we associate 
the force on a body with its rate-of-change of momentum. 
This notion of force is used consistently throughout this 
paper, but the reader is warned that it is common in the 
literature to instead label anything contributing to an 
acceleration as a force. 

It is useful to define a spin angular momentum vector 
S a associated with the center of mass frame. This is 
given by 

S a = - l -e a bcd {p h /m)S cd . (36) 

Since p a S a b = 0, S a b can be recovered from S a via 

Sab = £abcd(p c /m)S d . (37) 
The evolution equation (31b) now reduces to 

T)Q a 1 

-f- = (p a /m)(p b /m)S b - -e a bcd (p b /m)N cd . (38) 
as I 

Eq. (31a) may be used to replace the p b appearing here. 
The first term on the right-hand side is responsible for 
Thomas precession. It arises from the requirement that 
p a S a — remain preserved under time evolution. 

It should also be noted that the (rest) mass m does 
not necessarily remain constant under evolution. Using 
(31) and (33), it is straightforward to show that [12] 

^ = -z a F a + (p a /m)(p b /m)N ab . (39) 
as 

The right-hand side of this equation can be naturally 
transformed into a sum of total time derivatives and "in- 
duction terms" [12, 16]. 

We may now summarize how to find the motion of an 
extended body. Initial data must first be specified, con- 
sisting of the body's center of mass position z(S), linear 
momentum p a (s), and angular momentum S ab (s) — or, 
equivalently, S a (s) — at a chosen time s. The evolution 
of these quantities is then determined by the coupled first 
order system (31) and (35). Of course, these equations 
are of no use unless one has some knowledge of F a and 
N ab . At least in principle, the force and torque may 
be extracted from (27) and (32). As a practical matter, 
however, it is normally more useful to expand the force 
and torque using multipole approximations. The results 
of such a procedure are given through quadrupole order 
in (61) and (70) below. 

Lab frame centroid 

The problem of defining a center of mass worldline may 
be thought of as the specification of some notion of mass 
dipole moment that vanishes when evaluated about an 
appropriate origin. Our choices (33) and (34) accomplish 
this without making use of any special structures exter- 
nal to the body itself. The resulting worldline is known to 



possess a number of useful properties expected of some- 
thing declared a "center of mass" [12, 16, 18, 21, 32]. 

Various other possibilities exist if one allows the intro- 
duction of structures external to the body itself. Con- 
sider, for example, a spacetime endowed with a "pre- 
ferred" timelike vector field r a . This might be Killing, 
but more generally, it may be interpreted as defining a 
congruence of "lab frame observers." These observers 
have a 4- velocity field 

< = T a /V-nr b . (40) 

The "constant time" hypersurfaces £(s) used to define a 
body's momenta may now be generated by setting n a — 
nf^. Similarly, (34) may be replaced by 

nUz L )S ab (z L ,n L (z L ))=0. (41) 

The worldline formed by the set of all points zl satisfying 
this equation will be referred to as the body's "lab frame 
centroid." It is interpreted as the set of points about 
which the body's mass dipole moment vanishes according 
to the lab frame observers. 

It is instructive to compare the lab frame centroid of a 
given body with the actual center of mass defined by (33) 
and (34) above. In order to make this comparison, it will 
also be helpful to introduce a third, "hybrid" definition 
for a body's center, wherein we take n a = n£ as above, 
but now consider points zh satisfying 

P a (za, n L (z H ))S , afc(2 H , n L (z H )) = 0. (42) 

Momentum-velocity relations analogous to (35) are eas- 
ily derived for both the lab frame and hybrid centroids. 
In the hybrid case, there is no change other than placing 
subscripts "H" on all quantities appearing in that equa- 
tion. This means that (35) holds with the substitutions 
p a CM -+p a u = P a (z H ,n L ), Fq M ^F£ = F a (z R ,n L ), etc. 
Note that these quantities can differ from their counter- 
parts associated with the center of mass frame due to the 
different choice of hypersurface and origin. 

The momentum-velocity relation in the lab frame case 
may be obtained by directly differentiating (41). This 
yields 

(-npi)z£ = {-r h zt)pl - N?n - Sl b ztV c r b , (43) 

where p£ = p a (zL, ^l), etc. The apparent simplicity of 
this equation in comparison with (35) is due to the fact 
that it has not been manipulated into an explicit formula 
for i£. 

One may also define spin vectors associated with both 
the lab and hybrid frames. is naturally given by (36) 
with the usual replacements p^ M — > p^, etc. The form 
of the evolution equation (38) is clearly preserved in this 
case. S£ is, however, defined slightly differently. Let 

S£ = - 1 -e\ cd nlSi d . (44) 
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We then have 

= nlS h , b ztV c nl - -e a bcd nl(Nl d + 2p£ig). (45) 

The centroid equations (43) and (45) differ from their 
respective center of mass equivalents (35) and (38) 
mainly in the terms explicitly involving spin. At slow 
speeds (so i£ M — n£), (35) reduces to 

mCM^CM - PCM = (ScM X fcM)A"CM + ■ • ■ (46) 

where ". . . " includes terms nonlinear in the spin and 
terms independent of the spin (which are not necessarily 
negligible). The vector notation used in (46) denotes a 
decomposition with respect to an orthonormal triad or- 
thogonal to r a . The lab frame equivalent of this equation 
is obtained from (43), and reduces to 

m L z L - pl = -<Sl x a g + (47) 

where 

a\ = -n£V„n£ (48) 

is minus the acceleration of a lab frame observer. If r a 
were Killing, this would correspond to the "acceleration 
due to gravity" in a Newtonian limit. 

Eqs. (46) and (47) display a curious reciprocity. For 
a body moving in flat spacetime subject to an external 
force, there is a spin-dependent contribution to the hid- 
den momentum in the center of mass frame, while there 
is negligible spin-dependent hidden momentum in the lab 
frame (since, choosing r a Killing, the "acceleration due 
to gravity" is zero). On the other hand, for a body 
freely falling in a stationary, curved spacetime (i.e., a 
body moving under the influence of gravity), the situa- 
tion is reversed: there is negligible spin-dependent hid- 
den momentum in the center of mass frame (since the 
force F a approximately vanishes), but substantial spin- 
dependent hidden momentum in the lab frame. When 
gravity is treated as an external force, it contributes a 
spin-dependent hidden momentum of the correct form, 
but it does so in the "wrong" frame. 

The situation for spin evolution is in fact precisely 
analogous. At slow speeds, the timelike components of 
(38) and (45) may be written as 

n L , a ^^ = -(S CM -F CM )/m CM + ... (49) 
as 

and 

DS a 

n L , a ^r = S L ■ a g . (50) 

Again, only terms linear in the spin have been retained. 
The term on the right-hand-side of eq. (49) gives rise to 
the well-known Thomas precession. We see that a simi- 
lar reciprocity occurs with eqs. (49) and (50) as occurs 
with (46) and (47). In particular, gravity treated as an 



external force contributes a Thomas precession term, but 
does so in the "wrong" frame. 

Having compared the center of mass and lab frames 
at the level of time-derivatives of position and spin, we 
now directly analyze the (undifferentiated) positions zqm 
and zl as well as spins Sq M and ££. To do so, it is 
most convenient to first compare these quantities to their 
counterparts in the hybrid frame. We restrict attention 
to the case of a body moving in flat spacetime, so that the 
two different worldlines and associated quantities may be 
directly compared. For simplicity, we consider the case 
of no electromagnetic interactions, although the general 
external force density f a may be nonzero. 

We first compute the difference between zn and zl on 
a t = const, hypersurface in global inertial coordinates 
x* 1 = (t, x l ). We also choose the lab frame observers such 
that T a — d/dt. Now, a spatial shift in the origin clearly 
leaves the linear momentum unchanged. We therefore 
have P M = Pl = Pn- ^ i s useful to decompose this 
momentum such that 

P"=7¥(t< + F), (51) 

where is orthogonal to and 7 = (1 — ViV 1 )^ 1 / 2 . 
Using (lib) — which is equivalent to (24) and (25) in this 
case — it is straightforward to see that 

S£" = Sg" +2(zH-z h ) [lM P v] . (52) 

Contracting both sides with t v while using (41) and (42), 
this reduces to 

M(z h -z H ) = VxS lI . (53) 

Note that no slow-motion approximations have been 
made in this equation. Closely related results may be 
found in textbooks [25, 26] (in less general contexts). An 
equation essentially identical to (53) played a role in the 
post-Newtonian analysis of Landau-Lifshitz momentum 
flow during bobbing found in [8]. 

The difference between S£ and is also straightfor- 
ward to compute. Using (52), we find that 

S&=-Y{S£+T»ViSL). (54) 

The correction term here gives rise to the Thomas preces- 
sion of Sg. As previously noted, S£ does not experience 
Thomas precession in flat spacetime. 

We now compare the true center of mass zqm with 
zr. It is clear that zcm = for a free mass (i.e., one 
with external force density f a = 0) in flat spacetime, 
since the momentum is then independent of the choice 
of hypersurface. Somewhat more generally, suppose that 
f a =/= 0, but that the motion is nonrelativistic in the lab 
frame, i.e., \V\ 2 = ViV 1 <C 1. We further assume that 
T af a 0(|U||/|), as will be the case if the net exter- 
nal force F a is (nearly) orthogonal to the 4-momentum. 
Finally, suppose that the body's characteristic size R sat- 
isfies 

\f\R±/M < OQV\) « 1. (55) 
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This implies, roughly speaking, that the velocity does not 
change significantly during a light travel time across the 
body. It may then be shown that 



pcm = ^ + o(M 2 ), 



and 



(56) 



(57) 



S£ u = S£ + 0(\V\ 2 ). 
It follows immediately that 

Mzcm = Mz H + 0(|1?| 2 ) (58) 
Combining (54) and (57), we obtain 



S£ M = S£ + T»V t Sl + 0(\Vf 



(59) 

This accounts for how Sq M can experience Thomas pre- 
cession while S£ does not. Combining (53) and (58), we 
obtain 



Mz CM = Mz L + S L xV + <3(|Vf ) 



(60) 



when both worldlines are compared on a single t = const, 
hypersurface. This accounts for how there can be hidden 
momentum associated with spin in the center of mass 
frame while there is none in the lab frame. Although our 
derivation has been restricted to Minkowski spacetimc, 
a comparison of (46) with (47) suggests that a similar 
result also holds at least in stationary curved spacetimes 
(where r a may be chosen to be Killing, so has the 
interpretation of a "gravitational acceleration"). 



the conditions necessary to ensure an accurate multipole 
expansion may be found in [15]. 

Up to quadrupole order, the gravitational force and 
torque are [12, 14] 



pa 

grav 



1 J bcdf \/ a R 
6 



bedf 



M ab — j cd fl a jy b] 

lv grav — g ° dL cdf i 



(61a) 
(61b) 



where J abcd is the quadrupole moment of T ab , defined by 
dS r 'T a ' b ' aK l[cl 



[14, 23] 

jabcd 



2e b]d]fh H a , 



fZ h C 



(62) 



The notation here indicates independent antisymmetriza- 
tions of the index pairs (a, b) and (c, d). As noted before, 
Oa = V a <7(a/, z) denotes the derivative of Synge's world 
function and w c (x') is the time evolution vector field for 
the hypersurfaces S(s). Similarly, a aa > = V a V a 'fj(x', z) 
and 



H a a =[-a\,]-\ 



(63) 



where the "—1" denotes the inverse. We assume that the 
body is sufficiently small that this inverse exists on each 
cross-section S n W. Finally, Q abcd (x', z) is defined by 



C. Evaluation of forces and torques 

Returning to the general development where spacetime 
may be curved and F a b ^ 0, recall that the force and 
torque appearing in (31) can be extracted from (27) and 
(32) using explicit expressions for the generalized Killing 
fields in terms of their Killing data. This results in inte- 
gral expressions that may be found in Eqs. (12.10) and 
(12.11) of [14] (in the case f a — 0). These integrals can 
be expanded in a multipole series, and we now present 
the results of this up to quadrupole order. The formu- 
las in this subsection are valid for an arbitrary choice 
of worldline Z and vector field n a . However, one would 
not expect the quadrupole approximation to be valid un- 
less Z is close to the center of mass worldline and n a 
close to the choice (33). The metric and electromagnetic 
field must also vary slowly over each cross-section of the 
body. Note that this last assumption implicitly implies 
that self-fields are small, i.e., that the object under con- 
sideration is a test body 7 . A more precise discussion of 



7 As previously mentioned at the beginning of this section, modi- 
fying (24) and (62) in a particular way allows large self-fields to 
be included very naturally [9, 22, 23]. In this case, the laws of 
motion presented here still hold if only the "regular" Detweiler- 
Whiting component of the self-field is sufficiently small. This is 
normally much smaller than the full self-field. 



d) bb" 



e abcd {x',z)= / du(cr aa a (c a »<rp,a 
Jo 



\(v"{u),z) 



, (64) 



where y"(u) is an affinely-parameterized geodesic satis- 
fying y"(0) = z and y"(l) — x' . In flat spacetime, 



Qabcd g a ^ c g d ^ b 



(65) 



Note also that J abcd = J^ cd = J ab ^ and J^ d = 0. 

Following [18], J abcd may be decomposed into three 
simpler components via 

jabcd = jabcd + n [ajb]cd + n [cjd]ab_ 3n [ajb][c n d]_ (gg-j 

The tensors introduced here satisfy J ab = J^ ab \ J abc = 

^ja[bc\ ^jabed ^j[ab]cd jjab[cd] ^-q^ 



n Tab „ rrabc „ rrabed „ rrabed n. 

n b J — ribJ — ribJ — ridJ — U, 



abed 



(67) 



j[abc] = j[abc]d = q (g 8 ) 

J ab , J~ abc , and J abcd are interpreted as the mass, mo- 
mentum, and stress quadrupoles, respectively. 

Note that if the spacetime metric is Ricci-flat, the addi- 
tion of appropriately symmetrized multiples of the metric 
to the full quadrupole moment can be used to show that 
the laws of motion are unchanged if J abcd is replaced with 
an "effective quadrupole moment" of the form 



(69) 
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and J^ c have all of the same index symmetries and 
orthogonality properties as their non-subscripted coun- 
terparts, but are additionally trace-free. 

Up to quadrupole order, the electromagnetic force and 
torque are given by [14] 

Fkm = lF\i h - \Q bc ^ a F bc + \q bcd V a V b F cd (70a) 

Nf u = 2Q c ^F b ^ c - ^Q^ a {2V c F b ^ d - V d F b \). (70b) 

Q ab = is the body's electromagnetic dipole moment 
tensor, and Q abc = Q a ^ its electromagnetic quadrupole 
moment tensor. Definitions for these quantities may be 
found in Eqs. (8.11)-(8.13) of [13]. They satisfy 

n a Q abc = Q [abc] = 0. (71) 

The tensors Q ab and Q abc may be decomposed as fol- 
lows. Let 

Q ab = 2n^ a d b ^ - e abcd n c ^ d , (72) 

where n a d a = n a fj, a — 0. We refer to d a and /i a as the 
electric and magnetic dipole moments, respectively. We 
also write 

Q abc = i(3Q a[ Vl - Q abc ), (73) 

where Q ab = Q^ ab \ Q abc = Q a ^ bc \ and n b Q ab = 
n b Q abc = Q\ abc] = 0. Note that only the trace-free part 
of Q ab contributes to the force and torque if V a F ab = 0. 

In summary, the center of mass worldline of an ex- 
tended body may be found by solving Eqs. (31) and (35). 
Adopting the assumptions outlined above, the gravita- 
tional and electromagnetic forces and torques appearing 
in these equations should be well- described by (61) and 
(70), respectively. Note that this system is not com- 
plete in that we do not have evolution equations for Q ab , 
Qabc ^ Qr j abed _ ^his is as expected since such equations 
should depend on the type of material under consider- 
ation. They cannot be deduced from stress-energy and 
charge- current conservation alone. 

As a simple application of these approximations, a for- 
mula for the hidden momentum of an extended charge 
moving in curved spacetime can be straightforwardly ob- 
tained by substituting (61) and (70) into (35). Although 
there is no obstacle to including additional terms, we 
shall display the result here only in the approximation 
of neglecting all terms of quadrupole and higher order 
as well as all terms nonlinear in the dipole moments and 
spin. One immediately obtains 

mz a = Pa + (qe acdf F b c n d S f /m - 2Q c [a F b]c )n b . (74) 

The content of this equation becomes more clear at slow 
speeds. Use of (72) and the usual lab frame decompo- 
sition of F ab into electric and magnetic fields (see eq. 
(105)) gives 

mz = p + S x (qE/m) — ft x E — dx B. (75) 

Note that we have already encountered the second term 
on the right side of this equation in (46) above. 



III. SIMPLE EXAMPLES OF BOBBING 
MOTION 

We are now in a position to illustrate the presence 
of bobbing effects in a number of simple systems. First, 
Sect. Ill A analyzes the motion of two spinning balls con- 
nected by an elastic tether ( "tetherballs" ) in flat space- 
time. Test bodies moving in the electromagnetic and 
gravitational fields of cylindrically symmetric sources are 
then considered in Sect. Ill B. For the tetherballs, we will 
find that the bobbing is of an entirely "kinematic" nature. 
For the test body moving in the field of a line charge, we 
will find that there is also a "dynamical" contribution to 
the bobbing. Finally, for the test body moving in the 
field of a line mass, we will find that there are important 
differences in the analysis, but that the final results are 
very similar to the electromagnetic case. 

A. Tetherballs 

Consider two solid, uncharged masses A and B (the 
"balls" ) connected by a thin elastic band T (the "tether" ) 
in flat spacetime. The worldtubes Wa and Wb associated 
with the two balls are assumed to be disjoint. We also 
assume that the thickness of the tether is very small com- 
pared to the spatial extent of A or B (in any direction) . 
The total stress-energy tensor of the system is conserved, 

V b {Tl b + T]f + T% b ) = 0. (76) 

Consequently, if we view A, B, and T as a single system, 
the net force and torque vanish. Its center of mass moves 
at constant velocity in global inertial coordinates and its 
total linear and angular momenta are both conserved. By 
(35), the center of mass velocity of the complete system 
is collinear with its linear momentum. 

Nevertheless, the center of mass positions za(s) and 
zb(s) of the individual balls — as computed from and 
Tg b following Sect. II B — display interesting behavior. 
To simplify the problem, suppose that the tether is con- 
nected to A and B at their respective centers of mass 8 . 
It then follows from (27) and (32) that the tether exerts 
no torque on either ball: 

N ab = N ab = q ( 7? ) 

One also sees that the forces on A and B are given by 
appropriate integrals of f a = —V b T^ b . Using (31a), (35), 



8 It should be noted that since A and B undergo non-uniform mo- 
tion, they cannot remain "rigid." The material point in the body 
at which the center of mass is located may therefore change with 
time. Nevertheless, for motion that is not particularly relativis- 
tic and balls that are not particularly deformable, such effects 
should be very small. If necessary, we can assume that the ball 
is equipped with a servo-mechanism that continually adjusts the 
location of the tether's attachment to the ball to lie at the center 
of mass. 
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and (38), we find that the laws of motion for each mass 



are 



pa 



Dp a _ 
~d7 ~ 

•a a a 7716 c rid 2 

mz = p - e bcd F p S jm 

^=p»(F b S b )/m\ 
as 



(78) 
(79) 
(80) 



Subscripts "A" or "B" should be appended to each vari- 
able here as appropriate. 

Despite the relative simplicity of these equations, solv- 
ing for the general motion of the balls would be extremely 
difficult in general. The forces F£ and Fg are unknown, 
and modeling them would require analyzing the dynam- 
ical degrees of freedom of the tether itself. 9 In order to 
simplify the problem, we proceed as follows. Consider 
the global inertial coordinates (t,x l ) — which we refer to 
as the "lab frame" — for which the center of mass of the 
total system remains fixed at x % — 0. Use the timelike 
Killing vector r = d/dt to define "lab frame momenta" 
P M = Pl as m Sect. II B. For ball A, we then have 



Tfd 3 x, 



(81) 



along with similar expressions for ball B and the tether T. 
The spatial components of the total lab frame momentum 
vanish in the chosen frame: P l A + P B + = 0. 

We now assume that the balls are identical in their 
physical properties. We also also make the following ad- 
ditional assumptions: 

1. In the lab frame, the motion of the balls and tether 
is non-relativistic and (55) holds. 

2. Each component of the lab frame momenta satis- 
fies 10 P£ « P£. 

3. The motion of the system is symmetric with respect 
to rotations by ir about the x 3 -axis. 

Defining £ = d/dx 3 , our assumptions imply that P A + 
P£ ~ and P£C M = P^- Hence, 

PIC* = PbC - o. 



(82) 



9 In Newtonian mechanics, it would be consistent to assume that 
the tether has negligible mass and is arbitrarily stiff, so that it 
can instantaneously propagate forces from one ball to the other. 
Such assumptions would allow one to ignore the dynamical de- 
grees of freedom of the tether. However, in special relativity, 
if the tether satisfies the dominant energy condition, its tension 
must remain less than its mass density. It is would therefore be 
inconsistent to assume a tether of negligible mass exerting arbi- 
trarily large forces. For similar reasons, the tether also cannot 
propagate forces faster than the speed of light. 
10 This condition will hold if the mass of the tether is small, and 
that no portion of it moves much faster than the balls. 



It follows from (56) that the vertical momenta associated 
with the center of mass frame satisfy a similar equation, 



^=^-0(\V 2 \). 



(83) 



To lowest order in the velocity of the balls, the x 3 - 
component of the momentum-velocity relation (79) re- 
duces (for either body, with subscripts A or B dropped) 
to 



mz- ( — (S x F) ■ C/ra. 



(84) 



Since by (83) the left-hand side vanishes (i.e., £ momen- 
tum is conserved), this equation shows that the center 
of mass must move in the x 3 direction (i.e., "bob") to 
compensate the hidden momentum —S x F/m. If both 
masses are revolving around one another, Eq. (80) and 
the assumption of non-relativistic motion imply that the 
orbital timescale will always be much shorter than the 
time over which the spin experiences any significant evo- 
lution. Therefore, S may be taken as approximately con- 
stant for each individual orbit. In this case it is straight- 
forward to see that the bobbing occurs in phase with 
the orbit. The symmetry of the system implies that this 
effect is identical for both balls, so that there is a net 
bobbing of the whole orbital plane, just as in numerical 
simulations of spinning black holes. 

In the above analysis the bobbing of tetherballs is 
seen to be a direct consequence of the existence of the 
kincmatical hidden momentum present in (84). As dis- 
cussed in section II B, this kinematical momentum may 
be viewed as arising from the difference between the lab 
frame centroid and the center of mass. Along those lines, 
we can give an alternative, lab-frame-based analysis of 
the bobbing effect. The analog of (84) in the lab frame 
is simply P ■ C, — mz^ ■ £ (as follows from (43)), i.e., there 
is no hidden momentum in the lab frame. Then by the 
vanishing of the lab frame £ momentum (equation (82)), 
the orbit of the lab frame centroid is confined to a plane. 
However, the center of mass differs from the lab frame 
centroid by equation (60), and therefore bobs in and out 
of the plane. This is the approach to understanding the 
bobbing that was taken in the introduction. Whether 
viewed as a consequence of the U S x F" hidden momen- 
tum or of the difference between the lab frame centroid 
and the center of mass, we refer to this sort of bobbing 
effect as "kinematical bobbing" since it does not depend 
upon the details of the dynamical forces involved. It will 
arise in every example we consider. 

If we relax the symmetry assumptions and allow the 
masses to be unequal, otherwise similar assumptions lead 
to a net bobbing effect, 



(m A ZA + m B z B ) ■ C 




x Pa 



C (85) 



If both masses are comparable, the amplitude of the sys- 
tem's oscillation is of order for optimally-aligned 
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spin vectors. This corresponds at least qualitatively to 
what has been observed in simulations of spinning black 
holes binaries [2, 8, 33]. 

It is particularly interesting to consider the special case 
where two spinning masses are pulled (almost) directly 
toward each other. Suppose, for concreteness, that the 
two bodies A and B are initially placed along the a^-axis. 
Also let both of their spins initially be antiparallel and 
directed along (or against) the x 2 -axis. As both masses 
are pulled together by the tether connecting them, (84) 
implies that they will develop a net velocity along the 
x 3 -axis. If they experience an inelastic collision, conser- 
vation of total momentum implies that the resulting com- 
posite will have zero spatial momentum. It will also have 
no external forces acting on it, no hidden momentum, 
and therefore no final velocity. If the forces pulling the 
two balls together increase with decreasing separation, 
the transverse velocity of the two balls would be seen to 
first increase monotonically as they come together, and 
then reduce very quickly to zero as they collide. This 
peculiar sequence of events has been observed in numeri- 
cal simulations of two spinning black holes allowed to fall 
directly towards one another and coalesce [33]. 

B. Line sources 

We now present two further simple examples of bob- 
bing for test bodies moving in static, cylindrically sym- 
metric electromagnetic or gravitational fields. 

1 . Electromagnetism 

Consider the field of an infinite static line charge in 
flat spacetime. Let r a denote the unit time translation 
associated with the rest frame of this charge, and C, a the 
translational Killing field pointing along its axis of sym- 
metry. It is clearly possible to choose a vector potential 
A a = —tpr a satisfying C T A a = Cc_A a = 0. 

We now analyze the motion of a charged test body 
placed in this potential. Since ( a A a — and V a C fc = 0, 
it follows from (29) that 

Vi = PaC (86) 

is conserved. For simplicity, we restrict attention to the 
case where = 0. 

Although PaC a = 0, the test body will, in general, have 
hidden momentum. It can therefore develop a finite cen- 
ter of mass velocity along the axis of symmetry. Ne- 
glecting quadrupole and higher moments while retaining 
terms linear in Q ab and S ab , we find from (74) that 

mC a z a = C a n b {qe ab cd F cf S d n f /m - 2Q c [a F b]c ). (87) 

For slow speeds this reduces to (see equation (75)) 

mk-C=[Sx(qE/m)-flxE}-C (88) 



The term proportional to spin in this equation is of ex- 
actly the same form as the one appearing in (84), and 
gives rise to analogous "kinematical" bobbing. In this 
case, however, there is an additional pi x E hidden mo- 
mentum, which gives rise to a "dynamical" bobbing con- 
tribution as well. 

It is interesting to compare this dynamical hidden mo- 
mentum with the ^-component of momentum carried by 
the electromagnetic field itself, 

Vf M = J TH^x. (89) 

The purely electric field of the line charge carries no mo- 
mentum, while the momentum of the test body's field 
is (by assumption) negligible. Therefore, the only con- 
tribution to "pEM j g f rom "cross-term" between the 
line charge and test-body fields. In the limit of slow mo- 
tion, this may be computed in the manner used below to 
obtain (131), yielding 

Vf M = -(ft x E) ■ (. (90) 

Thus, the electromagnetic field momentum is equal and 
opposite to the dynamical part of the test body's hidden 
momentum. As the test body orbits the line charge, the 
electromagnetic field momentum oscillates 180° out of 
phase with the dynamical part the hidden momentum of 
the body, and in phase with the dynamical contribution 
to its bobbing velocity. 

It might appear that there is a failure of conservation of 
total ^-momentum in this case, since the ^-momentum of 
the body 11 is conserved by itself while the ^-momentum 
of the electromagnetic field oscillates. This is an artifact 
of our treating the line charge as "fixed." There is, in 
fact, a Lorentz force density exerted on the line charge 
by the electric field that results from the motion of the 
orbiting magnetic dipolc, which is straightforwardly in- 
tegrated to give a total force of d/di(/Z x E) ■ £. If the 
line charge were of finite length and thickness — so that it 
would have finite mass and one could meaningfully talk 
about its momentum — this force would impart to it a 
momentum (jlx E) • £ that exactly compensates for (90), 
restoring conservation of total momentum. 

The hidden momentum of the line charge may be 
shown to be negligible in this case, so, in fact, the line 
charge itself experiences a slight bobbing effect. 12 Per- 
forming a mass-weighted average of the £ _v elocities of 



1 On account of the symmetries in this problem, "Pc is in fact equal 
to the ordinary integrated stress-energy momentum f T 03 d 3 a; to 
the order of approximation relevant here. 

2 The moving line charge will now produce a magnetic field, so 
that its self-field may now carry some momentum of its own. 
However, we can ensure that this momentum remains negligible 
by taking a limit as the mass of the line charge becomes very 
big (at fixed linear charge density). Note that the mass times 
velocity of the line charge remains fixed in this limit. 



14 



both objects in this system results in an equation analo- 
gous to (85). The mass times bobbing velocity of the line 
charge exactly cancels the dynamical contribution to the 
mass times bobbing velocity of the orbiting body, leav- 
ing only the kinematical "S x a" term. The analogous 
conclusion for the case of two orbiting charges will be 
obtained more explicitly in Sect. V below. 

2. Gravitation 

Very similar effects also arise for an uncharged, 
spinning test particle moving in the metric of static, 
cylindrically-symmetric mass described by the Levi- 
Civita metric 

d s 2 = ~{x l fP&t 2 + (xy^-^Kdx 1 ) 2 + (dx 3 ) 2 } 

+ a- 2 (x 1 ) 2 - 4 "(dx 2 ) 2 . (91) 

x 1 is the cylindrical radius and x 2 an angular coordinate 
ranging from to 2tt. The central mass is oriented along 
the x 3 -axis. p and a are free parameters, the former of 
which is interpreted as a mass per unit length [34] . 

Consider the motion of a spinning test body in this 
metric. The translational Killing field £ = d/dx 3 has 
nonvanishing gradient, so it is the quantity 

V c = PaC + \s ab V a ( b (92) 

that is conserved rather than p a ( a by itself. Neglect- 
ing effects nonlinear in the spin or involving quadrupole 
or higher multipole moments of the body's stress-energy 
tensor, it follows from (74) that p a = mz a . If we consider 
the case Vq = 0, the vertical velocity of an orbiting test 
mass is given by 

mz a Ca = -\e abcd W a C b n c S d . (93) 

A straightforward calculation shows that the only non- 
vanishing components of V a £{, are 

ViCs = -VsCi = 2p(2p - lXx 1 ) 8 " 2 - 4 "- 1 . (94) 

Combining this equation with (93) determines the ver- 
tical velocity of the orbiting body. The result is more 
suggestive if written in terms of the fictitious gravita- 
tional acceleration a b g introduced by (48). Defining this 
using the timelike Killing field d/dt and using (93) and 
(94), we obtain 

mz a C = -(l-2p)e abcd n a <; b S c a d g . (95) 

The slow speed limit of this equation is 

mk-C = -(l-2p)(Sxa s )-C, (96) 

which is similar to (84) and (88) above. Once again, the 
bobbing velocity is found to be proportional to a term 
with the form "(spin) x (acceleration)". 



There are, however, several key differences from the 
previous two cases. First, linear momentum was (nearly) 
conserved for both the tetherball and electromagnetic 
examples. This meant that the bobbing could be un- 
derstood as arising from a time-dependent hidden mo- 
mentum. In the gravitational case, however, the hidden 
momentum is negligible. This is, however, compensated 
by the fact that it is not the linear momentum that is 
conserved by itself, but rather the particular combina- 
tion (92) of the linear and angular momenta. Second, 
it should be noted that the notion of acceleration used 
in (95) is not the true acceleration of the body — which 
nearly vanishes — but rather the acceleration of a station- 
ary observer. Finally, we note that p <C 1 in the "Newto- 
nian limit." If p is neglected, the gravitational bobbing 
here can be seen to have the same absolute value but op- 
posite sign to the kinematical part of the bobbing found 
in the tetherball and electromagnetic examples. 

In the next section, we will show that the close analogy 
between the electromagnetic and gravitational results is 
not an artifact of this particular example. Rules will be 
derived that allow us to translate equations for the mo- 
tion of a test body in a stationary electromagnetic field 
into those for a test body in a weak, stationary gravi- 
tational field. We will thereby account for the sign dif- 
ference in the S x a g bobbing of the test mass in the 
gravitational case as arising from the same kinematical 
term as in the tetherball and electromagnetic examples 
above, together with a dynamical effect that is a factor 
of two larger and of opposite sign. As in the electromag- 
netic case, the net bobbing of the complete system should 
be given only by the kinematical term. This illustrated 
directly by Eq. (125) below for the case of two compact 
masses in orbit around one another. 



IV. ANALOGIES BETWEEN 
ELECTRODYNAMICS AND GENERAL 
RELATIVITY 

It has long been known that there are strong simi- 
larities between electric and gravitational phenomena. 
This is clear when comparing Newton's law of gravity 
to Coulomb's electrostatics. In that case, the only differ- 
ences between the two theories lie in an overall sign as 
well as the existence of both positive and negative elec- 
tric charge. Many aspects of magnetic phenomena are 
also analogous to gravitational effects found in general 
relativity. Although distinctions between Maxwell's and 
Einstein's theories are extremely important, their many 
similarities can still be a useful source of intuition. 

We now review how — for the case of a stationary, 
nearly-flat spacetime with sources that are slowly mov- 
ing and have negligible stresses — the Einstein field equa- 
tions can be written in a form analogous to Maxwell's 
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equations 13 . We shall then consider test body motion 
and show that the laws of motion governing an extended 
mass in such a spacetime can be generated by precise 
rules from the laws of motion describing an object mov- 
ing in an electromagnetic field in flat spacetime. 

A. Field equations 

Consider a spacetime whose metric g a b satisfies Ein- 
stein's equation with stress-energy source t a b- Suppose 
that g a b is nearly equal to the flat metric rj a b, so the per- 
turbation h a b — g a b — r} a b will approximately satisfy the 
linearized Einstein equation with source t a b- We choose 
the Lorenz gauge, 

d b h ab = 0, (97) 

where 

Kb = Kb - ^VabV cd Kd, (98) 

and d a represents the covariant derivative compatible 
with r] a b. Unless otherwise noted, we raise and lower in- 
dices with the flat metric. The linearized Einstein equa- 
tion then reads [11] 

d c d c h ab = -167rto6. (99) 

The linearized Bianchi identity implies that t a b is con- 
served with respect to r] ab , i.e., d a t ab = 0. 

Let r a be a time translation Killing field of rj a b with 
unit norm with respect to that metric. This fixes a "lab 
frame." We define a mass-energy current j a with respect 
to r a by 

Ja = -t ab T b . (100) 

It follows immediately that d a j a = 0. 

We now restrict attention to the case of stationary per- 
turbations so that C T rj a b — H T Kb = 0. We further as- 
sume that the components of t a b projected orthogonally 
to r a are all negligible compared with the components 
contributing to j a . This means that in the lab frame 
adapted to r a , stresses in the body are all small com- 
pared with the local energy and momentum densities. In 
global inertial coordinates {t,x l ) where r = d/dt, it fol- 
lows from (99) and our stationarity assumption that each 



component of h^ u satisfies Laplace's equation with source 
-T67rf M „. Since £y has been assumed to be negligible, 
we have Kj — 0. The trace-reversed metric perturbation 
therefore has the form 

Kb = A( a 71b)cT C - (101) 

for some potential A a - As in electromagnetism, A a may 
be used to define a skew-symmetric field strength 

F ab = 2d [a A b] . (102) 

It immediately follows that 

d [a F hc ] = 0. (103) 

This is clearly analogous to one of Maxwell's equations. 

The other Maxwell equation is nearly, but not exactly, 
reproduced in linearized general relativity. Using (99), 
one finds that 

d b T ab = 32n\j a + ^r a (r b j b )} (104) 

instead of (8). Despite this, the field strength may still 
be split up into electric and magnetic components that do 
have precise electromagnetic counterparts. In Maxwell's 
theory, one normally defines electric and magnetic fields 
by 

F ab = 2r [a E b] + VabcdT c B d (105) 

and r a E a = r a B a = 0. We have used r\ a bcd here to denote 
the volume element associated with the flat metric. In 
the gravitational case, it is useful to define gravitoelectric 
and gravitomagnetic fields £ a and B a by 

F ab = S(T [a £ b] + VabcdT C B d ) (106) 

and r a £ a = r a B a = 0. Note the relative factor of two 
difference between the definitions of £ a and B a in terms 
of F a b as compared with those for E a and B a in terms of 
F a b- Up to terms nonlinear in the metric perturbation, 
the equations satisfied by £ a and Bb are now identical 
to the flat-spacetime Maxwell equations for E a and B a 
with source j a . It should be emphasized that this cor- 
respondence holds only for stationary fields. Fundamen- 
tally time-dependent phenomena like electromagnetic in- 
duction do not appear to have close analogs in general 
relativity [37]. 



A number of different comparisons can be made between Ein- 
stein's and Maxwell's equations. We choose to adopt a restricted 
post-Minkowski approximation here (see, e.g., [f 1] for a similar 
treatment of the field equations) . Alternatively, equations asso- 
ciated with the post-Newtonian approximation of general rela- 
tivity may also be successfully compared with equations found 
in ordinary electrodynamics [35]. It is even possible to derive 
some exact results by comparing various properties of the "tidal 
tensors" R a bc d an d Vo-Pjie [37] (see also [36, 38] for additional 
analogies) . 



B. Comparing laws of motion 

We now turn our attention to computing the laws of 
motion satisfied by freely- falling (F a b = f a = 0) extended 
test bodies in weakly-curved spacetimes like those just 
described. These will be compared with the equations 
governing test body motion in flat spacetime electromag- 
netism. The relevant external fields F a b and F a b will both 
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be assumed stationary and source-free near each body's 
worldtube. 

As a result of this analysis, certain substitution rules 
will be found that transform electromagnetic equations 
into gravitational ones. It is important to note that the 
direction of this correspondence cannot be reversed in 
any natural way. This is due to the fact that the elec- 
tromagnetic laws of motion depend on all multipole mo- 
ments of J a as well as the first two moments of T ab (the 
mass and spin). By contrast, motion in the gravitational 
case depends only on the multipole moments associated 
with the stress-energy tensor. It follows that the consti- 
tutive equations (and initial data) required to complete 
our laws of motion must be fundamentally different in 
the electromagnetic and gravitational cases. Even at the 
monopole level, for example, it would be difficult to take 
a single parameter — the mass — and uniquely transform 
it into both a mass and charge. 

An additional subtlety arises at dipole order. The 
dipole moment of T ab is essentially S ab , while the dipole 
moment of J a is Q ab . Both of these quantities can af- 
fect the motion of a charge in flat spacetime. In general, 
however, they may have very different structures. The 
center of mass condition (34) strongly restricts the form 
of S ab , for example. Physically, there is no need for such 
a constraint on Q ab . Despite this, simple substitution 
rules can only be found if we do restrict to special cases 
satisfying p a Q ab = 0. This corresponds to considering 
objects with vanishing electric dipole moment, and is the 
only case that will be considered here. 

We also introduce several approximations and special- 
ize to globally inertial coordinates aligned with the sta- 
tionary Killing field r = d/dt. Assume that the mass 
to, charge g, center of mass velocity z\ and linear size R 
of the test body all scale linearly with a small parame- 
ter 14 e. Internal (e.g., rotational) velocities are allowed 
to be of order unity, so we have S v ~ \x v ~ 0(e 2 ) and 
J7^p ~ c7tf ~ 0(e 3 ). Time derivatives of dipole and 
quadrupole moments will be neglected. 

We are interested in calculating the derivative 

£ = i a d a (107) 
ds 

of the body's rest mass to, spin S^, and "kinetic momen- 
tum" mi 1 , to order e 3 . Only terms linear in the fields E a 
and B a or £ a and B a will be retained. This calculation 
extends results obtained in [7]. Note that the derivative 
operator d a associated with r) a b appears on the right side 
of (107), so in the gravitational case, d^/ds = d 2 z l /ds 2 
measures the deviation from straight line motion in the 
metric r\ a b rather than the deviation from geodesic mo- 
tion in the metric rj a i> + h a \,. 



Note that this does not correspond to a post-Newtonian (or 
"post-Coulombic" ) expansion. That would take, e.g., m ~ 0(e) 
but i ! ~ O^ 1 / 2 ). See the next section. 



We compute d 2 z % /ds 2 in both the electromagnetic and 
gravitational cases by directly differentiating (35) while 
using (31a) to evolve the momentum. The forces and 
torques appearing in the resulting equations are those 
discussed in Sect. II C. In the gravitational case, these 
can be simplified in the weak-field, slow motion approx- 
imation being considered here. First note that to linear 
order in the metric perturbation, the Riemann tensor has 
the form 

Rabcd — ^^[adb^cd + T[ c dd]^ab) 

- d\a{nb][cFd)f)T f ■ (108) 

Substituting this into the formula (61a) for the gravita- 
tional quadrupole force yields 

*£«v =~\ J bcdf d a d c F d[f r lg]b T° . (109) 
Similarly, the torque (61b) reduces to 

ATgtv = lr df[a [rc(2d d ^ f d f ^ d ) 

-{T^d c F df -7 ]df d c T b \^)], (110) 

It is useful to replace both instances of J abcd in these for- 
mulas by the "effective quadrupole moment" J^ d with 
expansion (69) . We similarly replace the electromagnetic 
quadrupole components Q ab and Q abc in (70a) and (70b) 
by their trace- free components Q^f anc ^ Qtf- For the 
stationary source-free fields considered here, these sub- 
stitutions have no effect on the laws of motion. 

In the electromagnetic case, we find that the accelera- 
tion d 2 z t /ds 2 of a charged body with vanishing electric 
dipole moment in Minkowski spacetime is given by 

+ ■ (B - k x E)] + z^OAE x (S - ^-S)} 1 

TO 

+ \ Q&VdjEk - l -e klm Q^d l d 3 B m + 0(e 4 ). (Ill) 

This result corrects 15 Eq. (45) of [7] and extends it to 
include quadrupole forces and terms quadratic in z. Sim- 
ilarly, we find 

A [m -(-£.£)] = 0(e 4 ). (112) 



15 The term z^d 3 [E x(jl- ^S)Y in (111) arises from the "hidden 
momentum," and was not known to the author of [7] in 1972. It 
did not appear in the first edition of Jackson's text on electrody- 
namics, although part of it, namely z^dj(E X ft)', is mentioned 
in later editions [5]. The full term may effectively be found in a 
text by Suttorp and de Groot [28], although it does not seem to 
be widely known. In addition, a calculational error appears to 
have been made in Eq. (45) of [7]. 
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Thus, m is not constant. A constant of motion is ob- 
tained to order e 3 if we subtract from m the usual ex- 
pression, —jj, ■ B, for magnetic dipole energy. See [12] 
and [10] for further discussion. Also note that we are fol- 
lowing the quite confusing standard convention of using 
<9 l (/2 ■ B) to denote 1 '' ^d l Bj. Since jl is defined only on 
the worldline of the body, it makes no sense to take a 
spatial derivative of the quantity fl ■ B, but there is no 
good alternative vector notation for ^d l Bj. 

Very similar calculations may be used to find the coor- 
dinate acceleration of a test mass in a stationary, weakly- 
curved spacetime. This yields 

^dV = _ m r( 1 + 2 M2^ +4 J. x B-2(z-£)'zT 
as 1 

-2d i [S- (B- J*x £)} - z j d 3 {£ x Sf 
~ \ji k F d l d 3 £ k + leumjW&djB" 1 + 0(e 4 ), (113) 

which extends Eq. (44) of [7] to include quadrupole forces 
and terms quadratic in z. In this case, m remains con- 
stant up to 0(e 3 ). 

The evolution of the spin in the electromagnetic and 
gravitational cases can be obtained from (31a), (38), and 
the multipole expansions for the force and torque dis- 
cussed in Sect. II C. One subtlety in the gravitational 
case is that we are writing evolution equations in terms 
of a fictitious flat background metric r] a i>. The form of 
such an equation for will therefore be different from 
one for Sn = g^vS" . Of course, the same comment could 
also be made for the acceleration. There, however, it is 
clearly most natural to consider i\ It is less obvious that 
£ M should be more interesting than 5 p (or vice versa). 
We therefore write evolution equations for orthonormal 
frame components of the spin. In the electromagnetic 
case, we choose the lab frame basis ej = d/dx l , which is 
orthonormal in the metric i]^. In the gravitational case, 
we instead use 

e? = e?-^K x el (114) 

This is orthonormal to linear order in hnv with respect 
to the full metric g^ v . 

A direct calculation now shows that the spin of a 
charged particle in flat spacetime evolves via 

i(^ef ) = [p x (B - 'z x E)] t + ^(E- §)% 
as m 

+eiflQ^d k E l - ^Q&idjBk + 0(e 4 ). (115) 



The term z J dj[E X (ft — q/mS)] is to be interpreted similarly. 
However, if we had not neglected the time dependence of ft and 
S, this term would be replaced by z^d^E X (ft — q/mS)]. The 
derivative, i^du, along the worldline of E X (ft — q/mS) is well 
defined, and this derivative would act on ft and S as well as on 
E. 



Similarly, the spin of a small mass in a weakly curved 
spacetime satisfies 

^(Vf) = -2[# x (6-«f)], - (£■ S)ki 

-eijiJildkS 1 + ^iidjBk + 0(e 4 ). (116) 

Note that the forms of the Thomas precession terms 
q/m{E ■ S)z*i and —(£ ■ S)zi differ from the more familiar 
ones given in, e.g. [5] and [7]. This is because we work 
here with a tetrad aligned with the lab frame rather than 
rest frame 4-velocity. Up to this difference, the equations 
here agree with and extend those of [7]. 

By inspection, the following recipe allows one to trans- 
form the electromagnetic equations (111) and (115) into 
the corresponding gravitational ones (113) and (116): 

1. Identify all terms explicitly involving Si. Replace 
qEi by —mSi, and then leave as-is. 

2. To all other terms, send m —¥ m, q —¥ m, //* — > 
S\ E, -> Si, B t -► 2B it z l -► 2z\ Q% ¥ -► 
S^f -> 2J^p. Multiply the overall result by "-1" 
(excluding the spin terms set aside in step one). 

3. In the expression for acceleration in the gravita- 
tional case, add the term 2mi l (i J Bj)- 

The first of these rules essentially states that all spin- 
dependent terms in the electromagnetic equations are 
preserved in the gravitational case (with the force qEi 
replaced by the analogous force —m£i). These are ex- 
actly the same terms argued to be kinematical in Sect. 
II B, and the nature of this substitution further supports 
that interpretation. Our second replacement rule is more 
familiar. It states that dynamical effects in electromag- 
netism tend to have the opposite sign in gravity. Velocity- 
dependent dynamical terms also tend to be twice as im- 
portant in gravity as they are in electromagnetism (hence 

Bi -> 2B t , z l -> 2z\ and Q^f' "> 2J#). 

In addition to this, there is one term in the gravita- 
tional acceleration, namely 2mz l {zj£ 3 ) 1 that has no elec- 
tromagnetic counterpart. The gravitational and electro- 
magnetic equations cannot be made to agree perfectly at 
all orders, and the presence of this term signals that we 
have carried the analogy about as far as it will go (at 
least in this direction). 

Several points of caution should be noted about our 
substitution rules. First, they provide a recipe only for 
transforming one particular set of differential equations 
into another. Solutions to the gravitational equations are 
not necessarily related in any simple way to those of the 
electromagnetic equations. The given rules also apply 
only to the center of mass acceleration and spin evolu- 
tion equations. Momentum- velocity relations, forces, and 
torques cannot be reliably transformed in the same way. 
Even the center of mass acceleration and spin evolution 
equations will not necessarily correspond as we have out- 
lined if a different gauge is chosen for h^ v or a different 
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tetrad is used to decompose S^. Finally, the correspon- 
dence rules apply only to the equations obtained at the 
stated order of approximation. 



V. ELECTROMAGNETIC BINARY SYSTEM 

We now move to our electromagnetic analog of a 
black hole binary, which consists of two charged mag- 
netic dipoles with spin orbiting around one another in 
flat spacetime. Working in a slow-motion regime anal- 
ogous to the post-Newtonian limit of general relativity, 
we first analyze the bobbing motion of this system, and 
then explore the flow of momentum between the field 
and the bodies. Finally, we use this information to spec- 
ulate about kicks. At each stage in this process, we dis- 
cuss what conclusions may be drawn for the gravitational 
case. As in previous sections, 3-vector notation will be 
used to denote the spatial components of 4-vectors in 
global inertial coordinates. 



A. Laws of motion 

For ease of comparison with the gravitational case, 
we adopt a slow-motion limit of electrodynamics that is 
closely analogous to the standard post-Newtonian limit 
of general relativity (as applied to compact objects). In- 
troducing a small dimensionless parameter A, we scale 
each body's 3-velocity, mass, charge, spin, and magnetic 
dipole moment such that 

z l ~ O(VX), q ~ m ~ 0(A), S v ~ /*„ ~ 0(A 2 ). (117) 

Attention will also be restricted to the case of no electric 
dipole moment, and all quadrupole and higher moments 
will be assumed 0(A 3 ) or smaller. We also assume that 
d/Zj/ds ~ 0(A 3 ) (the same order of magnitude of changes 
in spin). Accelerations that contribute to the bobbing 
effect first appear at 0(A 25 ). Since this is 1.5 powers of 
A higher than the Coulomb acceleration, we will call this 
order of approximation "1.5 post-Coulombic," or 1.5PC, 
in analogy with terminology used in the post-Newtonian 
limit of general relativity. 

We determine the acceleration of each body by de- 
manding that it move as a test-body in the other body's 
field, according to the formalism developed in this paper. 
Taking into account that the electric field of each body 
is O(A) and the magnetic field of each body is 0(A 15 ), 
equations (31a), (70a), (72), (75), and (105) combine to 
give at 1.5 post-Coulombic order that 



■ q{l+^\zf){E + ixB) 

+ d i \p-{B-kx E)} 
d 



for the acceleration of each body caused by the electric 
and magnetic fields of the other body. Note the simi- 
larity of this equation to (HI), which is valid only for 
stationary external fields. In fact, equation (118) could 
have been derived from (111) by noting that the field of a 
body in its instantaneous rest-frame is stationary to the 
relevant order, using (111), and transforming to the lab 
frame. Finally, equation (118) also follows from the work 
of [9] and [10], where self-fields were properly taken into 
account. 17 

The first line of equation (118) is the ordinary Lorentz 
force. The second line is the ordinary dipole force (which 
too arises from the Lorentz force density). The third 
line corrects the integrated Lorentz force density for the 
presence of hidden momentum of the kinematical (first 
term) and dynamical (second term) variety. 

We now plug in the actual E and B fields to equation 
(118) in order to explicitly determine the acceleration. 
Formally, we begin with the Lienard-Wiechert field (and 
its analog for a magnetic dipole) and expand in powers 
of the velocity. However, only terms linear in velocity 
will survive at 1.5PC, so that retardation effects do not 
enter. We also specialize here to the terms in the accel- 
eration proportional to dipole and spin. (These terms 
are entirely responsible for bobbing, since in the absence 
of spin and magnetic dipole moment the orbit would be 
confined to a plane.) At 1.5PC order, the contributions 
to the acceleration of body A proportional to a magnetic 
dipole moment or spin are 



niAa A 



ftS -^[3(?xf)(ft 
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r 3 

QAQB 



v x ^ B J 



+ r [3((/Za x v) ■ r)r + 3(f • v)(Pa x r) 



2p A x 
(119) 



(5a xv — 3(v ■ f)(SA x f 

mA?" J L 

where r = Za — %b ls the orbital separation, and f = r/r. 
We have also used v = z A — Zb to denote the relative 
velocity of the two charges. The first line of (119) corre- 
sponds to the first line of (118), i.e., the ordinary Lorentz 
force on body A due to the dipole field of body B. The 
second line of (119) corresponds to the second line and 
final term of the last line of (118), i.e., the dipole effects 
on body A due to the Lienard-Wiechert field of body B. 
The final line of (119) corresponds to the first term of 
the final line of (118), i.e., the kinematical spin effect. 

We now focus on bobbing by considering only the ac- 
celeration perpendicular to the "instantaneous orbital 
plane" spanned by r and v. Let 



r x v 



(120) 



-S x (qE/m) +flxE] 



(118) 



17 In fact, it can be seen from these references that the Abraham- 
Lorcntz-Dirac self-force contributes to d 2 z l /ds 2 at 1.5PC order. 
This is unlike the gravitational case, where radiation reaction 
effects first appear at 2.5PN order. However, the self-force is 
always directed along the orbital plane and therefore cannot give 
bobbing. We therefore ignore the self-force effects on the orbit. 
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Since only terms proportional to spin and dipole can con- 
tribute to a A , this quantity may be computed from equa- 
tion (119). We find 



j_ QAf^B 

m A a A = 



r 3 |f x v\ 
qAqBS A /m A 



r°\r x v 



' /JJMA [f(-2v 2 + 3(f • v) 2 ) - v(f ■ v)] 
f(v 2 - 3(r • v) 2 ) + 2v(f ■ v)}. 



(121) 



This gives the acceleration out of the plane on each in- 
dividual body, thus showing that there is a "dynamical" 
contribution (involving /1a and /2b) in addition to the 
kinematical effect (involving S A ). As in Sect. Ill A, we 
now consider the mass- weighted average m A a^ + m B a^ 
of the bobbing accelerations, which provides a measure 
of the net oscillation of the system's orbital plane. This 
quantity is easily read off of equation (121) by noting that 
v and f are odd under the exchange A ■<-» B. This "an- 
tisymmetrizes" the kinematical contribution but in fact 
eliminates the dynamical contribution, so that the net 
bobbing is given by 



m A a A + m B a B — 



QaQb I S^A _ Sb_\ 
3 |f x v I I m A nig I 

f (v 2 - 3{v-f) 2 ) + 2v(v-r)] 



(122) 



Thus there is in fact no net dynamical bobbing in the 
electromagnetic binary: any dipole effect pushing one 
body upwards is always compensated by a corresponding 
dipole effect pushing the other body downwards. This 
cancellation is analogous to that found in section IIIB; 
the dynamical hidden momentum bobbing effect on one 
body is compensated by the Lorentz force on the other 
body. Although this cancelation is somewhat mysterious 
from the laws of motion point of view, in the next section 
its necessity will be seen at the level of conservation of 
momentum. 

Recalling that E B = q B f/r 2 + 0(A 15 ), equation (122) 
may be recast back into the canonical "5 x F" form of 
the kinematical bobbing effect (present already in (118)). 
This yields an expression for the net bobbing, 



toa^a + m B a B = 



d_ 

ds' 



Sb 
m B 



x (q A E B ) 



r x v 
If x v\ 



(123) 



which is directly analogous to the net bobbing result ob- 
tained in Eq. (85) for the tetherballs system. 



sufficiently stationary that the analysis of section IV for 
the motion of body A is valid to 1.5PC and 1.5PN orders. 
As in the electromagnetic case, the terms in the acceler- 
ation of each body that are linear in velocity can depend 
only on the relative velocity of the bodies. Therefore, the 
substitution rules of Sect. IV B should hold for equation 
(118) restricted to terms linear in the velocity. However, 
no non-linearities in velocity occur proportional to spin 
or dipole, so we can use our substitution rules on the 
spin and dipole terms of (118) and plug in the appropri- 
ate gravitoelectric and gravitomagnetic fields to obtain 
the analog of (121), 

m A a A = 

-2(m A S B + m B S A ) 



r°\r x v 
m B S A , 



[f(-2v 2 + 3(f ■ v) 2 ) - v(r ■ v)] 
f(v 2 - 3(f ■ vf) + 2v(f ■ v)]. (124) 



Although all terms are now proportional to spin, we 
would still say that the first line is dynamical in origin 
while the second line is kinematical in origin. The net 
bobbing again arises only from the kinematical terms, 



m A a A + m B a B - 

m A S B - m B S A 
r 3 |f x v\ 



[f (v 2 - 3{v ■ f) 2 ) + 2v(v ■ f )] . (125) 



This result agrees with the corresponding equation (5.19) 
obtained in [8] using post-Newtonian methods (after dot- 
ting the latter equation with f x v/\f x v\ and fixing a 
typographical error involving a missing factor of 1/r 3 ). 
Reintroducing the gravitoelectric field £ B ~ m B f/r 2 , 
equation (125) may alternatively be written 



m A a A + m B a B = 




x {- , m A £ B ) 



r x v\ 



(126) 



Equations (85), (122), and (126) are of an essentially 
identical form and clearly demonstrate the identical 
(purely kinematical) nature the net bobbing displayed 
by slow-motion binaries interacting via elastic, electro- 
magnetic, and gravitational forces. In each case, the net 
bobbing is controlled by the identical form of "5 x i 7 " ', 
where F is the particular force holding the bodies in or- 
bit. 



B. Momentum flow 



Gravity 

The corresponding results for a gravitational binary 
may be obtained as follows. As already remarked, in the 
instantaneous rest frame of (say) body B, the fields are 



Although the net bobbing occurring in our electromag- 
netic binary system is purely kinematical in nature, we 
now demonstrate that there is a considerable flow of mo- 
mentum between the bodies and the electromagnetic field 
during each orbit. 
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In order to calculate the momentum contained in the 
electromagnetic field, we note that the stress-energy ten- 
sor of the electromagnetic field is quadratic in F ab . 
It naturally breaks up into terms quadratic in the field 
of body A, terms quadratic in the field of body B, and 
cross terms linear in the fields of both bodies A and 
B. By the Schwartz inequality, the self-field contribu- 
tions to the electromagnetic stress-energy tensor cannot 
be negligible compared with the cross terms (which we 
shall denote by T x b ). Unfortunately, however, inclusion 
of these self-field terms would take us beyond the test 
body approximation, and therefore outside the scope of 
the formalism developed in this paper. We will simply 
appeal to the work of [9] and [10], which demonstrates 
that, at the order of approximation to which we work, 
stress-energy components quadratic in the field of body 
A can naturally be viewed as "self-energy" terms that 
comprise part of the stress-energy of body A. Similarly, 
terms quadratic in the field of body B can be viewed 
as comprising part of its stress-energy. "Renormalizing" 
both objects' stress-energy tensors in this way, the only 
self-field effect remaining in the laws of motion at this or- 
der is the Abraham-Lorentz-Dirac self-force. This makes 
no contribution to bobbing or to momentum flow per- 
pendicular to the orbital plane. All self-field effects will 
therefore be ignored here (as they have been above). 

We therefore focus attention on the momentum, P' x , 
carried by the cross-term portion, T x b of the electromag- 
netic stress-energy tensor, 

P* = J T° x \ (127) 

In vector notation, we have 

P x = ^J(E A xB B +E B x B A )d 3 x. (128) 

To the accuracy required for the 1.5PC approximation 
(which is 0(A 3 ) for the momentum), we may replace E A 
by — V$a, where $a is the instantaneous Coulomb po- 
tential of body A. Making a similar replacement for E B 
and integrating by parts, we obtain 




(129) 

It is now useful to perform the expansion 

$ A (2) = $a(.?b) - (x-z B ) ■ E A {z B ) + . . . , (130) 

where z B denotes the instantaneous center of mass of 
body B. Plugging this expansion into the integral of 
(129), using the definition of the magnetic dipole mo- 
ment, and keeping only terms proportional to dipole, we 
obtain 

^dipole = _ /7b x e a (z b ) - ma x E B (z A ). (131) 



Since only momentum proportional to dipole can point 
out of the orbital plane, this equation (dotted into the 
normal) gives the entire electromagnetic momentum per- 
pendicular to the plane. We have thus found that the 
relevant electromagnetic field momentum is in fact equal 
and opposite to the sum of the "/2 x E" portions of the 
hidden momenta associated with the two charges (see 
equation (75)). 

The reason for this remarkable-looking cancellation 
can be understood as follows. Consider the case where 
both the 3-current J A of body A and the charge density 
p B of body B vanish, so, in particular, fl A and q B vanish. 
There will be no force on either body if they are both 
stationary so this system can exist in stationary equilib- 
rium. Nevertheless, despite the fact that the bodies and 
their electromagnetic fields are stationary, the above cal- 
culation shows that the cross term in the electromagnetic 
stress-energy tensor still carries a momentum 

P x = -mb x E A (z B ). (132) 
However, the total stress-energy 

rpab rpab i rpab i rpab , rpab , rpab /-i qq\ 

J tot — J A + J B + J A,EM +i B,EM + i x K 1 ^) 

remains conserved, where is the stress-energy of the 
matter that makes up body A and Ta^em is the stress- 
energy of its electromagnetic self-field (and likewise for 
B). In particular, we have 

d T t ° o ° t + 9^ o ° t =0. (134) 

The first term vanishes by stationarity. Multiplying by x, 
integrating over space, and then integrating the resulting 
expression by parts, we obtain 

PL = J Tlotd 3 x = 0. (135) 

Note that neither the purely electric self-field of body A 
nor the purely magnetic self-field of body B can make a 
contribution to Pot- Consequently, (135) can be made 
compatible with (132) only by assigning a hidden me- 
chanical momentum fl B x E A to body B that is equal 
and opposite to the electromagnetic field momentum. In- 
deed, this way of resolving the apparent conflict between 
(135) and (132) was used in [30] to argue for the exis- 
tence of hidden momentum in electromagnetic systems. 
Note that in this stationary example, one may construct 
a simple model of a magnetic dipole to understand at a 
microscopic level the way in which the hidden mechanical 
momentum is stored by the dipole [6]. 

Given that exact cancellation of the electromagnetic 
field momentum and hidden mechanical momentum must 
occur in the stationary case, it should not be surprising 
that for quasi-stationary systems, there remains a can- 
cellation of the electromagnetic field momentum propor- 
tional to the dipole and the non-kinematical portion of 
the hidden momentum of the bodies. 
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The above considerations yield a rather surprising pic- 
ture of what happens in the 1.5PC approximation when 
two charged, spinning bodies orbit each other with their 
spins and magnetic dipole moments lying in the orbital 
plane. During an orbit, there is a substantial transfer of 
momentum perpendicular to the orbital plane between 
the electromagnetic field and the bodies. Nevertheless, 
when the momentum is transferred from the electromag- 
netic field to the bodies, it does not produce any net bob- 
bing motion; all of the perpendicular momentum trans- 
ferred to the bodies from the fields goes into the hidden 
momentum of the bodies, and there is no net effect on 
the perpendicular motion of their centers of mass. Net 
bobbing arises solely from the kinematic "a x S" con- 
tribution to the hidden momentum — exactly as for the 
tetherballs of subsection III A. 

As previously noted in subsection IIIB, the gravita- 
tional case differs on a fundamental level from the elec- 
tromagnetic one. The true acceleration (or force) and 
the true hidden momentum of each body (nearly) vanish 
in pure gravity. Furthermore, although the total mo- 
mentum of the complete system is well-defined if the 
spacetime is asymptotically flat, there is, in general, no 
meaningful way of assigning one portion of this total to 
the gravitational field and another portion to the bod- 
ies. Field momentum need not be localizable or additive. 
Nevertheless, the analogy established in Sect. IV sug- 
gests that in the 1.5PN approximation, each body acts 
as though it possessed a hidden momentum S x a (where 
a = d 2 z/ds 2 ). This may be interpreted as arising from 
the usual kinematical contribution — S x a as well as a 
dynamical term 2S x a (corresponding to \i x E in the 
electromagnetic case). Furthermore, the system acts as 
though, during each orbit, the dynamical portion of the 
hidden momentum is exchanged with gravitational field 
momentum. 



C. Kicks 

We have shown in this paper that bobbing behavior is 
a ubiquitous feature of the motion of relativistic spinning 
bodies. Furthermore, we have shown explicitly that for 
the tetherball system, for the electromagnetic binary at 
1.5PC order, and for the gravitational binary at 1.5PN 
order, the net bobbing is of a purely kinematic nature. 

In the case of coalescing black holes, not only is bob- 
bing behavior seen, but, in some cases, after merger the 
coalesced system receives a large "velocity kick." In some 
cases, this kick appears to be a continuation of the net 
bobbing motion of the black holes just prior to merger. 
Thus, it appears natural to expect that the velocity kick 
is closely related to the bobbing behavior and that, in- 
deed, these could be two aspects of the same basic phe- 
nomenon. We now argue that this is not the case, and 
that the kicks and bobbing are essentially unrelated. 

In the case of a merger of two bodies, the resulting 
body should move inertially and, by (35), would be ex- 



pected to have negligible hidden momentum. Thus, in 
order have a velocity kick, it is necessary for the merged 
body to have a momentum kick as well. In order for this 
to be possible, it is necessary that momentum be carried 
off by another system. In the case of tetherballs, this 
is impossible, since there is no other system available. 
Thus, although the tetherball system bobs in a manner 
very similar to the black hole system, if the tetherballs 
were to merge, the merged system could not acquire any 
velocity kick. 

On the other hand, a velocity kick should be possible 
in the merger of two charged spinning bodies. We have 
seen in the previous subsection that even in the 1.5PC 
approximation, a significant amount of momentum can 
be stored in the electromagnetic field (see (131)). If a 
merger event were to occur and if some 18 of this field 
momentum could be "released" and allowed to radiate to 
infinity rather than be reabsorbed by the merged object, 
then the merged object would acquire a velocity kick. 

The 1.5PC approximation is not adequate to analyze 
kicks, because at this level of approximation, the field mo- 
mentum (131) can only be exchanged with the bodies and 
cannot be radiated to infinity. However, it is not difficult 
to imagine that in a more relativistic system containing 
charged spinning bodies, a field momentum comparable 
to (131) could be radiated to infinity if a merger event 
were to occur. Although many of the quantitative details 
would undoubtedly differ from the 1.5PC approximation, 
we would expect (net) bobbing to be governed primar- 
ily by spin and kicks to be governed primarily by mag- 
netic dipole moment. In particular, if the charged bod- 
ies have spin but no magnetic dipole moment (e.g., they 
are composed partly of uncharged rotating matter and 
partly of charged non-rotating matter) then we would 
expect significant bobbing but no significant kick. By 
contrast, if the bodies have magnetic dipole moment but 
no spin (e.g., they contain oppositely charged portions 
that counter-rotate), then we would expect to have neg- 
ligible net bobbing effects, but in suitable cases it might 
be possible for a merged object to acquire a large kick. 
Similarly, if in an electromagnetic system we were to re- 
verse the magnetic dipole moments but not the spins, 
then we would expect the net bobbing to be the same 
but the kick to reverse sign. 

In the gravitational case, both the kinematic net bob- 
bing and the perpendicular field momentum are deter- 
mined by a x S at 1.5PN order and cannot be varied 
independently. Thus, in the gravitational case, it should 



If we perform the integral in (128) over a small ball enclosing 
body B rather than over all of space and then perform the same 
manipulations that led to (131), we find that a surface term con- 
tribution of j/Ib x now arises from the integration by parts. 
We thereby find that the field momentum contained within the 
ball is — | /Ib X E A ■ This calculation shows that 2/3 of the 
momentum (131) stored in the electromagnetic field is located 
within — or very near — the bodies, and would not likely be re- 
leased in a merger event. 
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not be surprising if the perpendicular kick velocity is cor- 
related with the phase of net bobbing at merger, and if 
the maximum perpendicular kick velocity is the same or- 
der of magnitude as the maximum net bobbing velocity. 
Nevertheless, it is our view that bobbing and kicks are 
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